Non-adiabatic Rectification and Current Reversal in Electron Pumps 
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Pumping of electrons through nano-scale devices is one of the fascinating achievements in the field 
of nano-science with a wide range of applications. To optimize the performance of pumps, operating 
them at high frequencies is mandatory. We consider the influence of fast periodic driving on the 
average charge transferred through a quantum dot. We show that it is possible to reverse the average 
current by sweeping the driving frequency only. In connection with this, we observe a rectification 
of the average current for high frequencies. Since both effects are very robust, as corroborated by 
analytical results for harmonic driving, they offer a new way of controlling electron pumps. 
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The possibility to pump charges or spins through 
nanoscale devices despite the absence of a bias voltage 
[U |2], shows the intriguing potential of driven nano- 
devices. It has consequently led to an increasing inter- 
est in such systems over the past decades. This interest 
was partly driven by the prospect of achieving single- 
electron pumping and thus creating a unique nano-scale 
single-electron source. However, it was soon realized that 
electron pumps may also help to close the metrological 
triangle, because they provide a connection between fre- 
quency and current [21 13] . These developments were sub- 
stantially aided by the rapid experimental and techno- 
logical progress in controlling and fabricating nano-scale 
devices. In particular, it recently became possible to real- 
ize charge pumping in the GHz regime [SHE] , which leads 
to a significant increase of the pumped current. 

The adiabatic limit of charge pumping, where the driv- 
ing frequency £1 is much smaller than the typical charg- 
ing/ discharging rate T, is very well understood [9| . In this 
case, one can resort to well-established time-independent 
formalisms to describe the electron transport and many 
theoretical works have considered adiabatic pumping in 
various systems [5HT5], The opposite limit of very fast 
pumping has also been studied Q31 [TS] ; to large extent in 
the context of photon-assisted tunneling [TBHTS] . How- 
ever, the borderland between these limits is lacking a 
comparably systematic understanding. Here, one is faced 
with an inherently non-equilibrium problem, while due 
to the similar time-scales a perturbative description is 
not possible. To address this problem, numerical calcu- 
lations in the time-domain [THl HO] or methods based on 
Floquet theory have typically been used [2TH21] . Only re- 
cently, within the diagrammatic real-time transport the- 
ory, has a summation to all orders in Q in the limit of 
weak tunnel-coupling and moderate pumping frequencies 
been achieved [T3] , revealing interesting non-adiabatic ef- 



fects for spin and charge pumping. 

In this Letter we demonstrate the drastic consequences 
of non-adiabatic driving for the pumped charge per pe- 
riod in a generic system and arbitrary frequencies. In 
particular, we show that it is possible to reverse the av- 
erage current by sweeping the driving frequency. This 
reversal is shown to be a hallmark of the borderland be- 
tween very slow and fast driving regimes. It does not 
rely on interference effects [32] and appears to be robust 
with respect to different driving schemes. The reversal 
effect is associated with a non-adiabatic rectification of 
the pumped current for high frequencies. 

In our case the electron pump is realized using a quan- 
tum dot (QD), which is coupled via tunnel barriers to 
source (S) and drain (D) contacts. These are connected 
to larger electron reservoirs. The total Hamiltonian is 
H = -ffdot+^rcs+-Htun, where the first term describes the 
QD itself, -ffdot = s(t)c'c, the second term characterizes 
the attached contacts, H ICS = I] QSS . D I] fes £afc&i fe fr Qfe , 
and the third term accounts for the tunnel coupling be- 
tween QD and contacts, H tUD = J2 a k T a k{t)b ak c + h.c. 
Here, c< and V ak create an electron in the QD state 
and in the reservoir state a k, respectively. The trans- 
port through the QD is mainly characterized by the tun- 
neling rate T a (t) — 27r^ fc \T ak (t)\S(e — e ak ), which is 
given in terms of the time-dependent tunneling ampli- 
tudes T ak {t). In order to calculate the pumped charge 
per period, we use the framework of NEGFs [25j [26] in 
conjunction with an auxiliary mode expansion |27j . This 
method is very flexible and allows for treating arbitrary 
time-dependencies of the parameters entering the Hamil- 
tonian. 

Specifically, we use the following time dependence for 
the QD level and the couplings to the reservoirs, 

e{t) = e + £i cos(fii) , (la) 

r a (t) = r° ex P [ Va cos(Oi - <y] . (ib) 
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This choice reflects the experimental situation of mod- 
ulated voltages, which will in general lead to an expo- 
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FIG. 1: Schematic representation of the pumping device: The level e and the couplings Tl,r to left and right reservoir oscillate 
in time, inducing periodic charging and discharging, a) Time dependence of e (black line, left axis) and Fl,r (red and blue 
lines, right axis) according to Eq. for a phase delay of 5=37r/4. b) Pumped charge Q per period as a function of frequency 
fl and phase shift 8 obtained from an NEGF calculation. The parameters defining the driving and the reservoirs are as follows: 
£ =0,£i = 20r, r) a = 6, r£ iR = re~ Va /2, fi a = 0, and k B T = T/10. 



nential dependence of the tunnel couplings on this mod- 
ulation [3IH]. The most important aspect of Eqs. (JlJ is, 
however, the inclusion of the phase shifts S a allowing for 
an offset of the coupling oscillations with respect each 
other. As in earlier work [7J, we use Si, = S and 5r = 
in the following. In Fig. [T^i the time dependence of e and 
r Q is shown for a specific delay of S = 37r/4. 

The time-dependent driving ([!]) induces currents Jl 
and Jr from the left and the right reservoir, respectively. 
The net charge, pumped from the left to the right reser- 
voir within one period r = 2tt/£1, can be obtained by the 
integral 



Q 



dt'lMt')- j R (t')} 



(2) 



In numerical calculations, the respective equations are 
propagated until the charge per period converges. Fig- 
ure [T] shows the numerical results for the pumped charge 
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FIG. 2: Pumped charge per period Q for the same parameters 
as in Fig. [I] NEGF (black lines) and rate-equation (broken 
orange) results are indistinguishable, a) Current reversal for 
the phase shift 5 — — 37r/4. b) Current rectification for the 
driving frequency fl = T/10. Parameters for which the time 
evolution is shown in Fig. [3] are marked by black points. Re- 
sult form the harmonic model [Eq. ( 12 I with ro = Ti = T/20] 



are shown as indigo/long-dashed line 



Q as a function of frequency ft and phase shift 8. For very 
low frequencies, f2 <C T, one finds the expected behavior 
of Q in dependence on 5: For negative shifts the pumped 
charge is positive, while for positive shifts it is negative, 
which is known as peristaltic pumping [S]. In striking 
contrast, one observes for higher frequencies (O > 10 _2 r) 
that the net current always flows in one direction. This 
implies for a negative phase delay (<$«— 37r/4) that by 
sweeping the driving frequency one can change the sign 
of the pumped current per period or, in other words, 
reverse the direction of the average current. These ef- 
fects — rectification and current reversal — are the cen- 
tral result of this Letter. 

In the following, we analyze the pumping using a sim- 
ple rate-equation description of the electron transport, 
which is valid for E\ 3> T, k^T. As we will show, this 
description is sufficient to reveal the basic mechanisms 
behind both effects. The currents Jl,r and the dot oc- 
cupation n are given by the following equations [7] 



Jl,rW = r L>R (t) [/(£(*)) - n(t)} 
d t n(t) = J L {t) + J R (t) , 



(3a) 
(3b) 



with /(e) the Fermi distribution function describing the 
occupation in the reservoirs. The pumped charge ob- 
tained in this model agrees very well with the NEGF 
result as can be seen in Fig.[2j which shows the current re- 
versal in panel a {Q vs frequency for 8 — — 37r/4) and the 
rectification in panel b (Q vs phase shift for — T/IO). 

In order to get a better understanding of this surpris- 
ing behavior, it is instructive to examine the temporal 
evolution of the currents for slow and fast drivings, re- 
spectively. This is most conveniently done using a Fourier 
analysis. By means of the definition ^ and the equation 
of motion ^ the pumped charge per period reads 



Q 



dt'T^(t')[f(e(t'))-n(l/)} 



(4) 



It is given in terms of the difference of the tunnel cou- 
plings r'~\ which is defined (along with the correspond- 
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FIG. 3: Time evolution within one period for three different parameters sets 5) . Lower row: Currents through tunnel 
barriers Jh(t) and Jn(t) (red and blue lines). Upper row: Dot occupation n(t) and transferred charge J gK'(JlA) — Jn(t'))/2. 
The gray area indicates the time span when the dot is charged, i.e., e(t) < 0. Red/blue arrows mark the times when left/right 
couplings are maximal. 



ing sum T^+A) which will be used later) as follows 



= r L (t)±r R (t) = £r 



(±) imnt 



A) 



Here, = r Q L I m ( Vh )e imS * ± T° R I m ( VK ) with I m is 

the modified Bessel function of the first kind of order 
m. The last expression is the Fourier series of the 
tunncling-rate sum/diffcrcncc. Analogously, the occu- 



pation n(t) = J2 T 



TL m C 



imQt 



and f{e(t)) = £ m im e 



hnQl 



can be expanded in a Fourier series. For sufficiently low 
temperatures, k^T <§; £i, we can replace the Fermi func- 
tion by the step function f(e(t)) = Q(—e(t)) and get 
f m = 9 m with 6 =l/2, 6 m =(-l)( m+1 )/ 2 /m7r for odd 
m, and 9 m =0 for even m^O. 

Plugging these series into Eqs. ^ yields an algebraic 
equation for the Fourier coefficients of the occupation, 
which reads in matrix-vector notation 



[0-n] 



(G) 



with the time-derivative operator -Dfc m = k8 m k and the 



coupling matrix Gkn 



-.(+) 



The components of the 



vectors 0, n and \ are given by the Fourier coeffi- 
cients introduced above [23]. It is important to notice, 
that neither D nor G contain the frequency O. This al- 
lows for a straightforward expansion in powers of il, as 
we will show below. For the charge Q in Eq. (H ) one needs 
n, which can be easily obtained from Eq. ( 6l) yielding 



Q=^f ( >. [0-n] , 



n = 



0. 



(7a) 
(7b) 



Note that this expression only depends on given quanti- 
ties, which are either external parameters (like f ^ ^ or Q) 
or trivial matrices (like D) . With this formulation one can 
derive intuitive expressions for low- and high-frequency 
pumping. In order to invert the matrix in Eq. Q we 
split the matrix G/il = Go + Gi into a diagonal and an 
off-diagonal component [29] 

(Go)fcm = 3km,Gkm/Q = ^km^o/^ 7 (8 a ) 
(Gi)km = (1 — 8km)Gkm/Q , (8b) 

and use for Eq. ^f]p) the expansion 







OO 

fc=0 







»-l<br J 



D-iGo)" 



(9a) 



(9b) 



where D— iGo can be easily inverted since it is diago- 
nal. Equation (JoJ) has a very intuitive interpretation. 
The Fourier vector D0 describes alternating "5-kicks" 
at times fitj = (7+1/2)77. The expansion in terms of 
accounts for the response of the system to these kicks, 
which is mainly characterized by the ratio Tq/ 51. The 
first term Kg in the sum is diagonal and given as 



(K ) r 



m + iro/51 

m 2 +(r /si) 2 



(10) 



The Lorentzian decay in the index m accounts for the 
exponential charging or discharging of the quantum dot. 
It is interesting to consider the following limits 



r /fi > 1 : 
r /o < 1 : 



n 



\Q — 1 



r Vr 



1A 1 + i-^fA 2 



(lla) 
(lib) 



4 



with Dkk = -Dfcfc for fc 7^ and £>oo = ~iFo/^, which 
replaces D in order to enable the matrix inversion. 

In the adiabatic limit {CI/Tq — > 0) given by Eq. (Hi), 
all matrix elements of Ko are identical. This implies that 
© n oc D0, i.e., the "5-kicks" of the driving mentioned 
before also occur in the response of the system. In other 
words, the electron in- our outflow is much faster than 
the external period. Indeed this can be seen in the lower 
panel of Fig. [3^. Therefore, Q in Eq. Q is determined 
by the couplings at specific times T^(tj). For negative 
phase shifts — 7r<<5<0, as shown in Fig.|3ji, it is 1^ '<0 
for the charging at Clto=-ir/2 and r^>0 for the discharg- 
ing at Clti=3ir/2. The opposite applies for positive shifts 
0<(5<+7r. As mentioned before, the pumping is "peri- 
staltic" 9J. The simple relation of Q and the coupling 
differences T^ito) and T^(ti) explains that the maxi- 
mal charge Q is obtained for S — ±tt/2 and that it van- 
ishes for 6 = 0, ±7r. Because of the pre- factor in Eq. Q, 
Q becomes independent of Cl in the adiabatic case and 
the first non-adiabatic correction is proportional to Cl. 

On the other hand, in the fast driving limit (CI/Tq — > 
oo) given by Eq. (flip), we get — n oc © + const and 
the occupation n(ijis constant. Because of Eq. Q there 
is no transfer in this limit and Q oc Cl . In between 
these two extrema the time scale of the exponential de- 
cay is comparable to the period of the external driving. 
Thus the integral Q "considers" r(~)(i) over the whole 
period, not just at particular instants of time as in the 
adiabatic case. This explains why the net current flows in 
the same direction as long as the peak of the left coupling 
occurs during charging of the dot, which is fulfilled for 
\5\>n/2, cf. gray area in Fig.[3] For <5=-3vr/4 (Fig.^) 
the charging occurs in the 2nd half of the charging pe- 
riod, for <5=+37r/4 (Fig.[3j:) in the 1st half. Thus, in both 
cases the dot is charged from the left and discharged to 
right since the right coupling Tr is locked to the oscil- 
lating level e. This explains the observed rectification 
effect. By means of this interpretation one would expect 
a current reversal for phase delays ~ir<5<—ir/2. For 
these delays it is r^^(io)<0, relevant for small CI, but 
J t 1 dt'r( - )(t')>0, relevant for large Cl, and the charging 
occurs either form the right or the left. Correspondingly, 
it is r(-)(ii)>0 and J t * 2 dtT(-)(i')<0 and the dot is dis- 
charged to the opposite direction. Figure [l] indeed shows 
this behavior in the predicted range of phase delays S. 

Finally, to show that the current reversal is not spe- 
cific to our driving scheme, we turn to the case of purely 
harmonic driving, i.e., T(t) = T + cos(Clt— 6). More 
importantly, the basic mechanism of the current rever- 
sal can be understood analytically in this case. Har- 
monic driving at a frequency Cl is characterized by having 



three Fourier components IY = {Fie +1(5 , To, Tie - " 5 } and 
Tr = {ri,ro,ri}, where Tq > Ti guarantees positive 
couplings and 5 is the time shift of the left and right 
coupling. If the calculation in Eq. (|9| is restricted to Ko 
only three Fourier components of the step function are 
needed, which are = {—1/ir, 1/2, —1/tt}. Using these 
expressions in Eqs. (|7|~(|To|) one gets 



Q = 



r 2 +ci- 



[Cl (l-cos<5) +T sin 6} 



(12) 



This simple expression contains all the basic features for 
periodic pumping including the current reversal. More- 
over, it allows us to analyze the respective regimes in 
detail. For Cl <C r (adiabatic limit) one obtains from 
Eq. (12 1 Q = (ri/r ) sin 5. As discussed earlier, the sign 
of Q depends on the order of the "door openings" . Op- 
timal transfer is attained for Ti = To. In the opposite 
limit, Cl > r , one finds Q = (Fi/0)(1- cos 6). Most in- 
terestingly and in contrast to the adiabatic case, in this 
limit the sign of Q is independent of the phase shift 5, 
which is the non-adiabatic rectification effect. Conse- 
quently, for S < one gets negative Q in the adiabatic 
and positive Q in the non-adiabatic limit: the average 
current can be reversed by tuning the driving frequency. 
These findings are confirmed by Fig.[2j where we compare 
the harmonic driving to the scenario considered initially 
[Eqs. ([T])]. Qualitatively, the behavior of Q is quite simi- 
lar in both cases, which underlines the robustness of the 
discussed effects. 

In summary, we have studied the influence of non- 
adiabatic driving on the charge pumping through a quan- 
tum dot in the Coulomb blockade regime. Our numerical 
calculations, based on a NEGF method, showed that the 
average pumped current can be reversed by sweeping the 
driving frequency. The origin of this effect was found 
to be the qualitatively different response to slow and fast 
driving, rendering the difference of the left and right tun- 
neling rates matter only at specific instants of time (adia- 
batic case) or during a time-interval (non-adiabatic case) . 
By means of a description with rate equations, we derived 
for the case of harmonic driving an analytical expression 
for the transferred charge per cycle, which confirms our 
analysis. Furthermore this shows that the observed ef- 
fects are generic and quite robust with respect to the spe- 
cific form of the external driving. Therefore they could be 
useful for realizing frequency filters or frequency-selected 
switches. 
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